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Abstract 

In the past decades, quantitative study of different disciplines such as system 
sciences, physics, ecological sciences, engineering, economics and biological 
sciences, have been driven by new modeling known as stochastic dynamical systems. 
This paper aims at studying these important dynamical systems in the framework of 
G-Brownian motion and G-expectation. It is demonstrated that, under the contractive 
condition, the weakened linear growth condition and the non-Lipschitz condition, a 
neutral stochastic functional differential equation in the G-frame has at most one 
solution. Holder's inequality, Gronwall's inequality, the Burkholder-Davis-Gundy (in 
short BDG) inequalities, Bihari's inequality and the Picard approximation scheme are 
used to establish the uniqueness-and-existence theorem. In addition, the stability in 
mean square is developed for the above mentioned stochastic dynamical systems in 
the G-frame. 
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functional differential equations 


1 Introduction 

Stochastic differential equations (SDEs) have been used profitably in a variety of fields in¬ 
cluding population dynamics, engineering, environments, physics and medicine. They are 
used to describe the transport of cosmic rays in space. For the simulation of the transport 
of energetic charged particles, SDEs solve several important equations such as Parkers 
transport equation and the Fokker-Planck equation [1], SDEs are also utilized in the field 
of finance as they are optimal to modern finance theory and have been broadly employed 
to model the behavior of key variables; the variables include the asset returns, asset prices, 
instantaneous short-term interest rate and their volatility. Biologists use these equations to 
model the achievement of stochastic changes in reproduction on population processes [2, 
3]. SDEs can be utilized to describe the percolation of a fluid through absorbent structures 
and water catchment [4]. There is a huge literature on the applications of SDEs in several 
disciplines of engineering such as computer engineering [5, 6], random vibrations [7, 8], 
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mechanical engineering [9-11], stability theory [12] and wave processes [13]. These non¬ 
linear equations in the framework of G-Brownian motion were studied by Peng [14, 15], 
Gao [16] and Faizullah [17]. The study of mild solutions for stochastic evolution equations 
in the G-framework was given by Gu, Ren and Sakthivel [18]. The pt\\ moment stability of 
solutions to impulsive SDEs in the framework of G-Brownian motion was established by 
Ren, Jia and Sakthivel [19]. In the G-framework, stochastic functional differential equa¬ 
tions were introduced by Ren, Bi and Sakthivel [20] and then generalized by Faizullah [21, 
22]. He also studied the p -moment estimates for the solutions to these equations [23, 24]. 
By virtue of linear growth and Lipschitz conditions, the existence theory for the solutions 
to neutral stochastic functional differential equations in the G-framework (G-NSFDEs) 
was given by Faizullah [25]. These equations not only depend on the present and past data 
but also depend on the rate of change of the past data [26, 27]. Subject to the non-linear 
growth and non-Lipschitz conditions, we do not know whether these equations admit so¬ 
lutions or not, if the solutions are unique or not and if they admit solutions whether they 
are mean square stable or not. The present article will fill the mentioned gape. Let 0 < to < 
t < T < oo. Suppose k : [to, T] x£C([-r,t 0 ];K rf ) — > R 1 *, X : [to, T] x£C([-r,t 0 ];M rf ) —> R d , 
p : [to, T] x 5C([-r, f 0 ];R d ) — *■ R rf and Q :£C([-r,t 0 ];M rf ) —> R J , are Borel measurable. 
We consider the following G-NSFDE: 

d[Y(t) - Q(t, Y t )} = <c(t, Y t ) dt + X(t, Y t ) d(B,B)(t) + /x(t, Y t ) dB(t), (1.1) 


where Y t - {Y{t + 9) : -r < 0 < 0, r > 0} indicates the collection of continuous bounded 
functions A : [—t, to] —*■ R rf equipped with norm ||A|| = sup_ r<6)<0 | A(0)|, which is a 
-8C([-r,f 0 ]; Revalued stochastic process and T(t) is the value of stochastic process at 
time t. All the coefficients K,X,p e M|([-r, r];R*) and {{B,B)(t),t > 0} is the quadratic 
variation process of G-Brownian motion {B(t),t > 0}. We denote by L 2 the space of all 
J^-adapted process T(t),0 < t < T, such that ||T || L 2 = sup^ r<{<7 . |T(t)| < oo. The initial 
data of equation (1.1) is JVmeasurable and a BC{[-r, to]; M rf )-valued random variable 


Y to = S = {m---r<0< : 0}, 


( 1 . 2 ) 


such that f to]; M^). In the integral form the above equation is expressed as 

Y{t) - Q{t,Y t ) = t;(0) - Q(to,Y to ) + f K {s,Y s )ds+ [ X(s,Y s )d(B,B)(s) 

Jtn Jtn 


■ / iAs, 


Y s )dB(s). 


An Revalued stochastic processes Y(t), t e [-r, T], satisfying the following features, is 
known as the solution of equation (1.1) with initial data (1.2). 

(i) The coefficients /r(t, Y t ) e £ 2 ([ o, T];R a ') and X(t, Y t ),/j.(t, Y t ) e £ 2 ([to, TJjR 0 '). 

(ii) Y(t ) is J^-adapted and continuous for all t e [to, T], 

(iii) To = f and, for each t € [to, T], 

d[Y(t) - Q(T f )] = K(t, Y t ) dt + X(t, Y t ) d(B,B)(t) + /At, Y t ) dB(t) q.s. 
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By a unique solution Y(t ) of G-NSFDE (1.1), we mean that it is equivalent to any other 
solution Z(t). In other words, we have to prove that 

e\ sup \Y(u) -Z(m)| 2 1 = 0. 

All through this article, the non-uniform Lipschitz condition, weakened linear growth 
condition and contractive condition are considered. These conditions are, respectively, 
given as follows. 

(Hi): Let t e [to, T], For every U, x £ -BC([-r, 0];!^) 

| K(t, U)-K(t,x) I 2 + I Ht, U) - k(t, x) I 2 + Wit, U ) - q(t, x)| 2 < T(| U-x | 2 ), (1-3) 


where the function T(- 
for w > 0 and 


/ 

J o+ 


dw 

TW 


= oo. 


R + is non-decreasing and concave with T(0) = 0, T (w) > 0 

(1.4) 


Since T is concave and Y(0) = 0, for all w > 0, 


Y(w) < a. + fiw, (1.5) 

where a and /S are positive constants. 

(Hu): For every t e [to, T] and ic(t, 0), \{t, 0), [i{t, 0) e L 2 , 

| K{t, 0) | 2 + | \{t, 0)| 2 + W{t, 0)| 2 < I<, (1.6) 

where K is a positive constant. 

(Hiii): For every [/,/€ £C([-r, 0];R J ) and t e [to, T], 

|Q(t, Y/)-Q(t, X )| 2 <« 0 |T/-X| 2 , 

, (1-7) 

|QM)| 2 <a 0 , 

where 0 < cuo < (g • 

Just for simplicity we consider d = 1 throughout the paper. 

2 Preliminaries 

This section is devoted to some basic literature, which will be helpful in our forthcoming 
research work. Firstly, we give three important inequalities. They are called Gronwall’s 
inequality, Bihari’s inequality and Flolder s inequality respectively [26]. 

Lemma 2.1 Let the function fit) be real and continuous on [a,b ]. Let K > 0 and (pit ) > 0, 
t e [a, b\. If fit) <K+ (fis)fis) ds,for every a <t < b, then 

fit)<I<ei- v(s)ds , 

for every a <t <b. 

The following lemma has been borrowed from [26, 27]. 



Faizullah et a\. Advances in Difference Equations (2017) 2017:350 


Page 4 of 14 


Lemma 2.2 Let the function T(-): R + —> K + be concave, continuous and non-decreasing 
satisfying T (0) = 0 and Y(v) > 0 forv > 0. Assume that<p(t) > 0, for all 0 < to < t < T < oo, 
satisfies 

4>{t)<K+ f <p(s)Y(fi>(s))ds, 

Jtn 

£ + . The following features hold. 

(ii) IfK > 0, we define U(t ) = ^ ds,for t e [to, T], then 

0(f) < U~ l (^U{C) + Jt ip{s)ds^, 
where Uis the inverse function ofU. 

Lemma 2.3 If for any q,r>\,- + -=\ and f,(p eL 2 then <pip e L l and 


where K is a positive real number and ip : [to, T] 
(i) IfK = 0, then 4>(t) = 0, t e [to, T], 

rt i 


f <P<P<(f 101 9 )Uf \<P\ T 

Ja \J a / \J a 


Next we specify some basic definitions and results of the G-expectation and G-Brownian 
motion [24, 28-31]. Let £2 be a nonempty basic space. Assume the space of linear real¬ 
valued functions defined on £2 is denoted by 'H. 

Definition 2.4 A functional E : U —> R is named a G-expectation if the following char¬ 
acteristics hold: 

(i) E[Y] < E[Z] whenever Y <Z, where Y,Z e H. 

(ii) E[y ] = y, where y is any real constant. 

(iii) E[0Z] = 9E[Z], for any 9 > 0. 

(iv) E[Y + Z\ <E[Y] +E[Z], where Y,ZeU. 

In the case when E holds only the characteristics (i) and (ii), then it is known as a non¬ 
linear expectation. For each w e £2 define the canonical process by B t (o>) = co t , t > 0. The 
filtration generated by the canonical process { B t , t > 0} is defined by T, = cr{B s , 0 < s < t], 
T - {.Ff}f>o- Let £2 be the space of all ^-valued continuous paths {w t ,t > 0} that start 
from 0 and B the canonical process. Also, suppose that associated with the distance given 
below, £2 is a metric space, 


p(w l , w 2 ) = — ( max I w] - w] I A1\ 

) Z^2'Vt 6 [0.A:] 1 f f| I 


Let Ci.Lip(R* rxaf ) denote the set of bounded Lipschitz functions on R Arxaf . Fix T > 0 and set 

L% (£ 2 r ) = {<P(B h ,B t2 ,...,B tk ) :k>l,h,t 2 ,...,t k € [0, T],f e C b . Up {^ kxd )}, 

£°,(£2 f ) C L°,(£2r) for t < T and L° ip (Q) = The completion of £?,(£2) under 

the Banach norm £[| • \ p ]p , p > 1, is denoted by L P G {i 2), where L P G {i2 t ) C L p g (Q.t) C L p G {i 2) 
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for 0 < t < T < oo. Suppose tzt = {to, h,..., t^}, 0 < to < fi < • • • < tx < oo is a partition of 
[0, T], Set p> 1, then M G °(0, T ) indicates a collection of the following type of processes: 


N~1 

>1t(w) = ^2 %i(w)I [ti ,t M ]{t), (2.1) 

i =0 

where f,- € L P G (Ei ti ), i = 0,1 ,...,N - 1. Furthermore, the completion of Mg°(0, T) with the 
norm given below is indicated by M P G { 0, T),p> 1, 

M = {f £ [l'hl p ]ds 

Definition 2.5 The canonical process {£(t)} f >o under the G-expectation E defined on 
L° ip (Q.) and satisfying the following properties is called a G-Brownian motion: 

(1) 5(0) = 0. 

(2) The increment B t ^ - B t , for any t, k > 0, is G-normally distributed and independent 
from B tv B G ,...,B tn , for n € N and 0 < t\ < t 2 < ■ ■ ■ < t n < t. 

In the G-framework the Ito integral /()/) and the corresponding quadratic variation pro¬ 
cess {(B)j}t>o are, respectively, defined by 

T N -1 

I(ri) = / r) u dB u = V] 5/(5f i+1 - B ti ), 

;=o 

(B) t =B 1 t -2 ( B u dB u . 

Jo 

The capacity c(-) associated to a weakly compact collection of probability measures V is 
given by 

c(D)= sup P(.D), DeB{Q), 

PeV 

where B(Q) is the Borel a-algebra of Q. D is known as polar set if it has zero capacity, 
that is, c(D) = 0 and a characteristic holds quasi-surely in short (q.s.) if it is valid outside 
a polar set. For each Y € L°(Qt)> the G-expectation E is given by E[Y] = swp PeV Ep[Y], 
where for each PeP, Ep[Y] exists. For more details of the above definition, we refer the 
reader to [28]. 

3 Some important lemmas 

Firstly, we define the Picard iterations sequence as follows. Let, for t e [to, T], Y°{t) = t, (0) 
and Y p = t; for each k = 1,2,..., then 

Y k (t) = Q(t,Yf)-Q(t,Y? o ) + k( 0)+ [ k(s, Y k ~ v ) ds + f k(s, T^ 1 ) d(B,B)(s) 

J to J to 

+ [ n{s,Y k ~ l )dB{s), 

Jto 



t e [to, T], 


(3.1) 
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In the following lemma we show that Y k (t) e r];R /c ). Then we prove another 

important lemma. These lemmas will be used in the existence-and-uniqueness theorem. 


Lemma 3.1 Suppose that assumptions Hj, Hu and Hi„ are satisfied. Then, for all k> 1, 


sup E\ Y k (t)\ 2 < a, 

-x <t<T 


ts- 10 boi* it * \ 

wherea = p^e 1 - 15 ^ 1 ' to) ,K 2 = [l + 15a 0 + 10ba*{T-t 0 )]Ef \ 2 +K h I< x = 5[l + 3ar 0 ]£|C| 2 + 
60ao + 10 K(T - to)a*{K + a), a* = q-i + a 2 + 0:3, aq, a 2 and a 3 are positive constants. 

Proof It is obvious that Y°{-) e M|([-r, T]; W k ). Using the basic inequality \c\ + c 2 + C 3 + 
c 4 + C5I 2 < 5|ci| 2 + 5\c 2 \ 2 + 5 | c 3 | 2 + 51C41 2 + 5|cs| 2 , equation (3.1) follows; 


Y 1 


it) | 2 < 5|?(0)| 2 + 5|Q(t, Y k ) - Q(t 0 , ^o)| 2 + 5 |/ ^( 5 ’ Y s _1 ) ds 
If k{s,Y k - l )d(B,B)(s) +5 f p(s, Y k_1 ) dB(s) 

I j to Jto 


+ 5 

By taking G-expectation on both sides, using Lemma 2.3 and the BDG inequalities [16] 
we obtain 

£|r*(t)| 2 < 5£|c(0)| 2 + 5E|Q(f, Y k ) - Q(t 0 , Y k o ) | 2 + 5 ai E f |/c(s, Y k ^)\ 2 ds 

Jto 

if \k(s,Y k ^)\ 2 ds + 5a 3 f \ii(s,Y k -')\ 2 

Jto Jto 


(*t ~ ft 

+ 5 ol 2 E I | A. (s, Yf 1 ) | “ ds + 5a 3 / \ij.(s, Y k ~ l ds 

Jto Jt 0 

5£|C(0)| 2 

+ 5E\Q(t, Y k ) - Q(t, 0) + Q(to, 0) - Q(t 0 , Y*) + Q(t,0) - Q(t 0 ,0)| 2 
+ 10o'i E f (|/c(s, Y* -1 )-/r(s,0)| 2 + |k(s,0)| 2 ) ds 

Jto 

+ 10 a 2 E f (|A(s,Y s A: - 1 )-A(s,0)| 2 +|k(s,0)| 2 )J(s) 

Jt 0 

+ 100-3 f (|/z(s, Y* -1 ) -/u(s,0)| 2 + |/r(s,0)| 2 ) d(s) 

Jt 0 

5£|^(0)| 2 + 15E|Q(f, Y k ) - Q(t, 0)| 2 + 15£|Q(t 0 ,0) - Q(t 0 , Y*) | 2 
+ 30£|Q(t,0)| 2 + 30£|Q(t 0 ,0)| 2 +10a 1 £ f |k(s,0)| 2 Js 

j to 

+ 10o , i£ f |k(s, Y* _1 ) - k{s, 0)| 2 ds 

Jt 0 

+ 10 a 2 E f \k{s,0)\ 2 ds + l0a 2 E f |k(s, Y^ 1 ) - A(s,0)| 2 d{s) 

Jto Jto 

+ 10o?3 f |/x(s, 0) |~ d(s) + 100:3 f \p,(s, Y* _1 ) - /u.(s, 0)| 2 d{s). 

Jto Jto 
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By using assumptions H n and Hin, we have 


E\Y k (t)\ < 5£|t:(0)| +15a 0 £|if | +15a 0 £|itr + 60«o 

+ lOcei K{T -1 0 ) + 10 a 2 I<(T- t 0 ) + 10a 3 K{T - t 0 ) 

+ 10«i e[ r(\Y s k - 1 \ 2 )ds + 10a 2 E f r(\Y s k ~ l \ 2 )ds 
J to J to 


'to 


+100*3 [ rdrf 1 ! 2 )^) 

Jt 0 

= 5£|f (0)| 2 + 60ffo + 10 I<(T - t 0 )(ai + 0:2 + a 3 ) 

+ 15ao£|lf| 2 + 15Q'o£|l^| 2 + 10(ai+a2 + a 3 )£ f Y(| if- 1 ] 2 ) ds 

't 0 

< 5£|f (0)|“ + 60«o + 10 I<(T - to)(ai + 0:2 + a 3 ) + 10a(ai + a 2 + a,^)(T - to) 

C t 

+ 15a 0 £| Y k \ 2 + 15« 0 £| 4T | 2 + 10j6(q!i + a 2 + « 3 )£ / \Y k ~ k \ 2 ds 

Jt 0 

= 7<i + 15a 0 £ \Y k \ 2 + 10Z7(q'i +a 2 + a 3 )E f | Y k ~ l \ 2 ds, 

J to 


where 7<i = 5£|f(0)| 2 + 60«o + 10 K{T - to)(ai + ot 2 + a 3 ) + 10u(ai + o ?2 + a 3 )(£ - to) + 
15ao£|(| 2 . We note that 


sup|lf| 2 < sup \Y k (q)\ < |(| 2 + sup \Y k {q)\\ 

to<s<t -T<q<t to<q<t 


from which follows 


sup £|l^)| 2 <£|C| 2 +/<i + 15a 0 sup E\Y k (q)\ 2 

-r <q<t ~Y<q<t 


+ 10/1 («i + a 2 + a 3 )£ f sup \Y k 1 (q)\f ds. 

Jt 0 -*<q<t 


For any j > 1, we observe that 


max£|yf 1 |“ < £|7| 2 + max£l Y A (s)| , 
1 <k<j 1 1 1 <k<j 1 1 


and we obtain 


max sup E\Y k {q)\ 2 < £|(| 2 +/<i + 15ao sup E\Y k {q)\ 2 

l<k<j- T <q<t -T<q<t 


+ 10d(«i + «2 + « 3 )E [ r£|(| 2 +max sup £|£ A '(^)| 2 1 ds 

Jto *- l<k<j -T<q<t j 


to 

: I< 2 + 15a 0 sup E\Y k (q)\ 2 

-r <q<t 


10P(ai + a 2 + a 3 )£ / max sup £ I Y k (q) I 2 ds, 

Jto -T<q<t 
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where Ki_ = £|£| 2 + K\ + 10/3(ai + a 2 + 0 : 3 )(T - to)E\k\ 2 . We have 


max sup E\Y k (q)\~ < 

l<k<j -z<q<t 


I <2 10 / 8(«1 + 0?2 + 0 : 3 ) 


1 - 15«o 


1 - 15«o 


/.t 

/ max sup £|Yf( 7 ) ds. 

Jt 0 l<k<j- T <q<t 


Consequently, the Gronwall inequality gives 


max sup £|l^(f )| 2 < a, 

1—^—7 -r <q<t 


10^(g 1 +«2 + <«3) (j- t j 

where a = 1 _ 1 g aQ e 1_15 “° , but / is arbitrary, so 

sup £|T*(f )| 2 <a. 

—x<t<T 

The proof is complete. □ 

Lemma 3.2 Let hypotheses Hi, H;i and H^ hold. Then for every k, d > 1, a constant j3 > 0 
exists such that 

e\ sup \Y k+d {s)-Y k {s)\ 2 ]</3 f T (e\ sup | Y k+d ~\q) - Y k ~\q)f\) ds 

< y(t-to), 


where y = /SY(4C) and f = 4(a '^“^ Q ' 3) . 

Proof By the basic inequality |ci + C 2 + C 3 + cff < 41 Ci | 2 + 41C 2 1 2 + 41C 3 1 2 + 41C 4 1 2 , equation 
(3.1) yields 


\Y k+d {t)-Y k {t)\ 2 <4\Q(t,Y k+d )-Q(t,Y k ) 2 + 4 f [, K(s,Y k+d ~ 

Jto 

+ 4 

f [A.(s, Y k+d - 1 ) - k(s, if' 1 )] d(B,B)(s) 
\Jto 

+ 4 

f t [ / i(s,Y k ^- i )- / 4s,Y k - l )]dB(s) . 

\Jt 0 


Next, on both sides, we take sub-expectations. Then we apply the Jensen inequality 
E(w(z)) < w(E{z)), BDG inequalities [16] and hypotheses Hj, Hu, H„i to find 

e\ sup |T* +d (s)-T*(s)| 2 l 

L-r <s<t -I 

<4a 0 £[ sup \Y k + d (q)-Y k {q)'f\ 

I—r <q<t 

+ 4cq f T (e\ sup \Y k+d -\q)-Y k -\q)f\)ds 

Jto ' '--E<q<s 

+ 4 a 2 f T (e\ sup | Y k+d -\q)-Y k -\q)\ 2 ])ds 
+ 4a 3 f r(4 sup | Y k+d ~ 1 {q)-Y k ~ 1 (q)\ 2 V\ds 

Jto V l -r<q<s 
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<4a 0 £[ sup \Y k+d (q)-Y k (q)\ 2 ] 

• — T<q<t -I 

+ 4(a 1 + a 2 + a 3 ) f Y^eT sup | Y k+d ^(q) - Y k ^(q)\ 2 V\ ds, 
e\ sup \Y k+d (s)-Y k (s)\ 2 ]<p [ t(e\ sup \Y k+d - 1 {q)-Y k ^{q)'\'\\ds, 

' — T<s<t -I Jto ' ' — r<q<s 

where p = . Finally, we use Lemma 3.1 and get 

E\ sup \Y k+ \s)-Y k {s)\ 2 }<pr{4C){t-t 0 ) = y{t-t 0 ), 

L-r <s<t J 


where y = pY(4C). The proof is complete. □ 

4 Existence-and-uniqueness results for G-NSFDEs 

In this section first we construct a key lemma. We set 


fiit) = y{t-t 0 ), te[t 0 ,T], 


(4.1) 


Let us define a recursive function as follows. For every l,d> 1, 

fk+iit) = p f Y (f k [s)) ds, 

Jt 0 

<l>k,d{t)=E\ sup | Y k+d {q)-Y k {q)f\. 

*--x <q<s 

Select Ti e [to, T] such that 


(4.2) 


pr(y(t-t 0 ))<y, (4.3) 

for all t e [to, T\\. 

Lemma 4.1 Let hypotheses Hi, Hu and Hm hold. Then, for all d > 1 and any k>l,a posi¬ 
tive Ti e [to, T] exists such that 


0 < 4>kAt) < fk(t) < fk-i(t) < ■ < fiif), 


(4.4) 


for all t e [to, T], 


Proof Mathematical induction is used to prove the inequality (4.4). We use Lemma 3.2 
and the definition of the function (/>(•) to obtain 


fi 4 {t) = E\ sup \Y ud {q) -Y 1 (q)\ 2 ] < y(t-t 0 ) = fi{t), 

'--z<q<s 


M)=E\ sup \Y 2+d (q)-Y 2 (q)\ 2 ]<P f t(e\ sup | Y l+d (q) - Yfq)\ 2 ]) ds 

'—T<q<s J Jto ' '—T<q<s J ' 


<P 


f r (Ms)) 

Jto 


ds = f 2 [t). 
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By (4.3), we obtain 


(t) = p[ T(0i {sj)ds= ( /3r(y{t-to))ds<y{t-t 0 ) = 4> 1 {t). 
Jtn Jtn 


Hence for every t e [to, T{\, we derive that <p2,d{t) < flit) < </>i(t). Assume that (4.4) is 
satisfied for some k > 1. Then we need to verify that the inequality (4.4) holds for k + 1. 
We proceed as follows: 


f k+u (t) = E\ sup \Y k+d+ \q)-Y k+ 1 (q)\ 2 ] 

'—r <q<s -I 

</l^ T (e[ sup \Y k+d (q)-Y k (q)\ 2 Jjds 

= pf T (4> k4 (s))ds 
Jt 0 

<p f T((l) k {s))ds 
Jto 

— 0/r+l(^)- 


Also 


4 >k + i(t) = p f T ( 4 > k (s))ds<p f T( 0 ir-i(s)) ds = 4 > k (s). 

Jtf) J tc\ 


Hence for all t e [to, T{[, 4> k +i,d{t) < (f> k +i{t) < 4> k {s), that is, Lemma 4.1 is true for k + 1. The 
proof is completed. □ 


Theorem 4.2 Suppose hypotheses H,, Hjj and Hm are valid. Then there exists at most one 
solution of the G-NSFDE (1.1) having initial data (1.2). 


Proof Firstly, we derive the uniqueness of solution. Let the G-NSFDE (1.1) with initial 
condition (1.2) admit two solutions T(t) and Z(t). Then we get 


\Y(t)-Z(t)\<\Q(t,Y t )-Q(t,Z t )\+ [ \K(s,Y s )-K(s,Z s )\ds 

Jto 

+ [ \Hs,Y s )-Hs,Z s )\d(B,B)(s)+ \ \yi{s,Y s )-ii{s,Z s )\dB(s). 

J tft J to 


On both sides, we take G-expectations and use the basic inequality (ci + C 2 + C 3 + C 4) 2 < 
4(c^ + c\ + c\ + cf). Then applying the Holder inequality and the BDG inequalities we get 


£|T(t)-Z(t)| 2 <4£|Q(t,T t ) 


Q(t,Z £ )| 2 + 4«! [ E\k(s, Y s ) - /c(s,Z s )| 5 
Jto 


ds 
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Using assumptions Hj, H,, and Hi,, we have 


l\ sup |r(^)-Z(^)| 2 l < 4a 0 £r sup j Y(q) - Z(^)| 2 1 

'—x <q<t -I ■—x <q<t 

+ 4(o?i + a 2 + a 3 ) [ r(£[ sup \Y(q) - Z(q)\ 2 ~\) ds, 

Jto V L -T <q<t J/ 


and it follows that 

E 


\ sup |r(^)-z(^)| 2 l 

'—x <q<t -I 


2 "I 4(«1 + a 2 + a 3 ) 


1 - 4ao 


f t(e\ sup \Y{q)-Z{q)\^ 1 \)ds. 

Jt 0 V L -r <q<t J/ 


Consequently, Lemma 2.2 yields£[sup T<(?<{ | Y(q)-Z(q) | 2 ] = 0, t e [to, T], The uniqueness 
proof is completed. To show existence we note that <p k (t) is continuous on t e [to, 7)] and 
decreasing on t e [to, 7i] for k > 1. We now use the dominated convergence theorem to 
define the function <p(t) as follows: 




lim <p k {t) 

k^-oo 


lim/S f T((f>i ( _i(s)) ds = ft f T (6(s))ds, 
k ^°° Jto Jto 


to < t < '/'i. 


So, 


<m < 0 ( 0 ) + p 



ds. 


Hence for every to < t < T\, Lemma 2.2 yields <p(t) = 0. For all t e [to, Ti], from Lemma 4.1 
it follows that <pk,d(s) < <pk(s) — > 0 as k — > oo, which gives E\ Y k+d (t) - Y k {t) \ 2 —> 0 as k —> 
oo. Then from the completeness of L 2 and assumptions H,, Hj,, H,;, follows that, for all 
t e [to, Ti], 


Q(t, if) -> Q(t, T f ), Jc(t, Y t k ) -> K(t, T f ), L(t, if) -* A(t, T t ), 
q.(t,if) —> /i{t, Y t ) ini 2 ask —> oo. 

Hence, for all t e [to, Ti], 

lim Y k (t) = lim Q(t, if) - Q(t, ff (0) + lim f k(s, if -1 ) ds 

k—>oo l —>■ oo k^-oo J 

+ lim f X(s,Y k ^) d{B,B)(s) + lim f ii(s,Y k ~ l ) dB(s), 
Jto k ^°° Jto 


that is, 


T(t) = Q(t,T f )-Q(t,C) + ?( 0 ) + 


[‘k(s, 

Jto 


Y s ) ds 


f X(s,Y s )d(B,B)(s) + f i4s,Y s )dB(s). 

J to J t(\ 


Hence the G-NSFDE (1.1) having initial data (1.2) admits a unique solution Y(t) on t e 
[to, Ti]. By iteration, we see that equation (1.1) admits at most one solution on t e [to, T]. 
The proof is completed. □ 
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5 Mean square stability 

In this section, we study the mean square stability for stochastic dynamical system (1.1). 
The following definition is borrowed from [26, 27]. 

Definition 5.1 Let Y(t) and Z(t) be any two solutions of the G-NSFDE (1.1) having the 
respective initial conditions f and £ belong to M 2 ([- r, 0] : K ( ). A solution Y ( t ) of equation 
(1.1) having initial data (1.2) is known to be mean square stable if for every £ > 0 a 5(e) > 0 
exists so that £|f -£| 2 < 5(e) implies E\Y(t) - Z(t) | 2 < e for every t > 0. 

Theorem 5.2 Suppose hypotheses H, and H„ are satisfied. Let equation (1.1) admit two 
solutions Y(t ) and Z{t) with initial data £ and f, respectively. Let t e [0, T], Ififior all e > 0 
a 5(e) > 0 exists such thatEf -£| 2 < 5(e), then 

E\Z(t)-Y(t)\ 2 <e. 

Proof Let system (1.1) admit two solutions Y(t ) and Z(t). Then, for any t e [0, T], it follows 
that 


Y(t) = $(0)-Q(t 0 fi) + Q(t,Y t )+ [ K(s,Y s )ds+ [ X(s,Y s ) d{B,B)(s) 

J t() J to 


+ f p (s, Y s )dB(s), 

J to 


Z{i) = m-Q{tofi) + Q{t,Z t )+ f i<(s,Z s )ds+ f k(s,Z,)d{B,B)(s) 

J to j to 


+ f n(s,Z s )dB(s). 

JtQ 


Then 


Y(t) - Z(t) = m- m - Q(to, 0 + Q(to, ?) + Q(t, Y t ) - Q(t,Z t ) 


+ 


+ 


f [K(s,Y s )-K(s,Z s )]ds+ [ [Hs,Y s )-k( S ,Z s )]d{B,B)(s) 
Jto Jto 

/ [/r(s, Y s ) - p{s,Z s )\dB{s) q.s. 

Jt n 


We use the fundamental inequality (cy + Ci + C 3 + C 4 + C 5 + cfi 2 < 6 (c 2 + c\ + c\ + c\ + c\ + 
Cg) and the sub-expectation on both sides. Then using the Holder inequality and BDG 
inequalities [16] to obtain 


?[ sup \Z(r)-Y(r)f]<6Ef-l;\ 2 

l— x <r<t -I 

+ 6a 0 £r sup |Z(r) - Y(r) | 2 ] + 6a 0 Ef-^ | 2 

L -r <r<t - 1 

+ 6(cti + a 2 + (X3) f T (e\ sup |Z(r) - T(r)| 2 ~h ds. 

Jto ' L _ T < r < 5 J/ 
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It follows that 

E\ sup |Z(r)-r(r)| 2 l < 6(1 + ao) Z|f; ~ II 2 

— r<r<t J 1 — 6a 0 

Finally, using Lemma 2.2 we get 
E[\Z(t)-Y(t)\ 2 ]<e, 

for t € [0, T], The proof is complete. □ 

6 Conclusion 

Neutral stochastic functional differential equations (NSFDEs) play a key role in model¬ 
ing physical, technical, biological and economic dynamic systems such as predicting op¬ 
tion pricing and the growth of populations. In general, one cannot obtain the explicit so¬ 
lutions to NSFDEs. The study of properties and behavior of solutions to NSFDEs such 
as uniqueness, existence and stability require extensive observations. In this article, we 
have used some useful inequalities such as the Flolder’s inequality, Gronwall’s inequality, 
the Burkholder-Davis-Gundy (in short BDG) inequalities, Bihari’s inequality and the Pi¬ 
card approximation scheme to obtain the existence and uniqueness of a solution for neu¬ 
tral stochastic functional differential equation in the G-framework. Moreover, the mean 
square stability is developed for the above-mentioned stochastic differential equations. 
The G-Brownian motion theory is not based on a particular probability space and gener¬ 
alizes the classical Brownian motion theory in a non-trivial way. The methodology used 
in this article to estimate the existence, uniqueness and stability of solutions for NSFDE 
in the G-framework is attractive and applicable in numerous practical applications. For 
example, the above-mentioned theory is useful in distributed system control [32], eco¬ 
nomics [33] and biological as well as neural control systems [34]. This article will play a 
key role in providing a framework for future work in this direction such as the study of the 
p-moment estimates for NSFDEs driven by G-Brownian motion and existence theory for 
stochastic pantograph differential equations driven by G-Brownian motion. 
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